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VIA. ON CLASSICAL OPTICAL INTERFERENCE, COHERENCE AND 
FLUCTUATIONS 

Light from a real physical source is never, strictly speaking, monochromatic nor does it 
emanate from a single point in space. Both the amplitude and the phase of the wave field 
generated by a real source undergoes irregular fluctuations. Within a time interval of the 
order of the putative coherence time the amplitude and phase of the radiation remain 
relatively stable and, thus, it may be treated as monochromatic. For time intervals long 
compared to the coherence time, the effects of fluctuations become manifest. In particular, 
interference phenomena, the hallmark of wave behavior, are observed only when the 
interfering beams traverse paths that differ in length by an amount small compared to the 
coherence length — i.e. the product of the coherence time and the velocity of light. In 
this first section we review the classical analysis and interpretation of interference effects in 
terms of the notions of partial coherence . 

Complex Representation of Polychromatic Fields: 

Let us first carefully define the mathematical framework of the subject. We assume that the 
observable optical field may be represented by a Fourier integral transform 1 




[ VIA-1 ] 



Since of necessity E' (t) must be real 



For this assumption to be valid, it is sufficient that the function E{t) be square-integrable — i.e. 
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£(-co) = E*((o) [ VIA-2 ] 

and, consequently, Equation [ VIA-1 ] may be re-expressed in terms of an integral over 
positive frequencies - viz. 

E (r \t)= |{£(co)exp[-/co?] + c.c.} d(0 ■ [VIA-3] 



Further, if we write 

£(q)) = ^ a(co)exp[/^G))] , [VIA-4] 

where a(co) and <f)(co) are both real, then the observable optical field may be expressed in 
the form 

E (r \t) = J a(to) cos[/ [(^(co) -co t]}d(0 I VIA-5 ] 

o 

which is a strict generalization of the complex or phasor representation of real 
monochromatic fields. To complete the generalization, we follow Dennis Gabor 2 and 
introduce the so called complex analytic signal 

E(t)= J a(co) exp{* [<^(co)-co?]}dco = 2 J £(co) exp [- i co t] d(0 [VIA-6] 



which is to be associated with the real observable field - viz. E (r \i) = Re[E(t)\ _ As we 
shall see, this rather careful, or should we say pedantic, introduction of the analytic signal 



2 D. Gabor, /. Inst. Elec. Engrs. (London), 93, 2529 (1946). 
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greatly facilitates and validates the use of complex quantities in handling cycle-averaged 
properties ~ e.g. energy and intensity — of polychromatic optical fields. 3 In most 
applications of interest in optics, the spectral amplitudes will have appreciable values only 
within a spectral width Aco which is small compared to some mean frequency (co) . Thus, 
it often useful to express the analytic signal in terms of a temporal modulation of a mean 
monochromatic - viz. 

E(t)= a(t)exp{i[§(t)-{(0)t]} . [VIA-7] 
Using Equation [ VIA-6 ] we see that 

a (Oexp [i$(tj\ = 2 J E (co) exp[-z' (co- (co)) f] d(0 = J E (co) exp [-i(Ot] d(o[ VIA-8 ] 

-<a>) 

where £( ro ) = 2 e[(0+ (co)) : by assumption, £(co) will be appreciable only near co= (co) . 
The following list of transforms are included for reference: 

a (t) = ^[E ir \t)} 2 + [E {i \t)} 2 = jE(t)E*(t) = \E(t)\ [ VIA-9a ] 



3 The following identities involving improper functions are invaluable when treating analytic signals — see Beran 
and Parrent: 

^- JexpJ/ x(y'-y))dx= §{y'-y) 
J- |exp[T^(/-,)] dx= 5 ± (/-y) = i L 5(/-y)+l-^-j 
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§(t)= (co)?+tan 




,r E\t)-E{t) i 

[ l E*(t)-E(t)\ 



[ VIA-9b ] 



E [r \i) = a(t) cosfz [${t) - (to) t]} 
E [i \t) = a(t) sin{/[4)(?)-(0))?]}. 



[ VIA-9d ] 



[ VIA-9c ] 



Using Equations [ VIA-1 ] and [ VIA-6 ], we may easily establish a form of Parseval's 
formula 



We have thus far assumed that field is defined for all values of t. In practical terms, the 
field or, more likely, the observation of the field is defined only within some finite time 
interval -T <t <T . In all instances of any consequence, T is much larger than the 
physically significant times 2tc/Aco and 2ti/(q)) so that we may feel some confident in 
taking the idealization T— > oo . However, stationary 4 requires that the time averaged 
energy of the field, which is proportional to 





[VIA- 10] 








A random process is characterized as stationary when all of the moments of the random variables are 
independent of absolute time. 
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T 



lim — 

t^2T 



dt ' 



[VIA- 11 ] 



-T 



must be finite and, thus, a stationary field function is not square-integrablel The problem 
of analyzing such functions has been a much discussed issue in the mathematics literature. 5 
However, as applied physicists, we brush these mathematical niceties under the rug in what 
follows. We proceed by firmly asserting that optical fields of interest do, indeed, have 
Fourier transforms and chance the consequences of trucation errors. 

The cycled- averaged intensity or power spectral density is a crucial element in the analysis 
of coherence effects and, from Equation [ VIA-6 ], it is proportional 



Consideration of this expression leads us to a logical definition of the first-order 
temporal correlation function of the electric field as 




[VIA- 12] 




T 




[VIA- 13] 



-T 



So that Equation [ VIA- 12 ] becomes 



5 See, in particular, N. Wiener, Acta. Math.,55, 1 17 (1930). 
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|£(w)| 2 =^j ) (E(t+x)E*(t))exp[i(i>%] dx= ^ j j r(x) exp[* cox] dx [VIA- 14 



and the normalized power spectral density (distribution) is defined as 



r 



F((o)= |£(co)| | |£(co)| da 
it J 



V 



lr, 2 »-|-i 

, (E(t+x)E*(t)) exp[zcox] dx \ (\E(t)\ ) 
1% jL\ /_ 



1 



— I y (1) (x) exph'on] dx 
2% J 



[VIA- 15] 



where 



Y (1) 



(E(t^)E*(t)) 
{E(t)E?(t)) 



[VIA- 16] 



This normalized correlation function is usually called the field's complex degree of 
first-order temporal coherence . 6 

Young's Interference Experiment -- A Rudimentary Measurement of 
Temporal Coherence 

By considering in some detail the simplest kind of experiment in which coherence effects 
are important — i.e. Young's (or Michelson's) observation of interference fringes — we can 



6 It is easy to show that 



F[(a) = (l/7t) Re | y (1) (x) exi{ i m ] dx 
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identify the significance of a light source's degree of first-order temporal coherence and 
demonstrate how that coherence can be measured. Consider the following basic 
experimental configuration: 



Source 




Colli mati ng 

Lens 

Screen 1 

The time integrated/averaged intensity at a given point P is proportional to 




l(P)=2(jE^P,t)] 2 ) = (E(P,t)E\P,t)) 



I 



Observation 
Point 

Q 



Screen 2 



[VIA- 17] 



The real instantaneous field the point Q is derivable from the following analytic signal 

E(Q,t) = Q E(P x ,t- t l ) + C 2 E(P 2 ,t- t 2 ) 

= C 1 E(P 1 ,t-sJc) + C 2 E(P 2 ,t-s 2 /c) ' 



[VIA- 18] 



It follows that the time integrated/averaged intensity at the point Q is proportional to 
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/(G) = QC?{E (p l3 t- t)E*{P v t- 0)+QC* (E(P it t- t)E*(P 2 ,t- t 2 )) 

+C 2 C*{E{P 2 ,t- t 2 ) £*(7>,*- t)) + C 2 C* 2 {E{P 2 ,t- t 2 ) E*{P 2 ,t- t 2 )) ' 

If the field is stationary , we may with confidence shift the time origin in the various 
expressions so that 



[VIA- 19] 



[E{P n ,t-t n ) E*{P n ,t-t n ))=(E{P n ,t)E*{P n ,t)) 

= /(*.)= 12(0) 

(E{P n ,t- t n ) E?(P m , t- t m ))= (E (P n ,t+ t m - E*{P m ,t)) 



r (1) (t - 1) = r 

nm \ m n) t 



The normalized coherence function 



[ VIA-20a ] 



[ VIA-20b ] 



[ VIA-21 ] 



is often called the complex mutual degree of coherence . Finally, we write the 
general interference law for stationary optical fields as 



7 (e) = 7,(2) + 7 2 (e) + 2 JtTq) JIJq) y« 

where I T {Q) = \q( 7 (/>) and I 2 {Q) = \C 2 ( I (P 2 ) . 



5 2 *2 ^ 



V c ) 



[ VIA-21a ] 
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Thus the degree of first-order temporal coherence may be expressed directly in terms of 
well-defined experimentally observable quantities — viz. 



where l{Q) and I„(Q) are, respectively, proportional to the intensity measured at Q with 
both pinholes open or with only the nth pinhole open. The visibility of the fringes 



Models of Chaotic Light Sources 

A Collision-Broadened Source - An example of homogeneous- 
broadening: 

Perhaps, the simplest model of a chaotic source is that of a collision interrupted 
oscillator. In this model one pictures an ensemble of atomic radiators each of which 
emits a field of constant amplitude, oscillating at a fixed frequency (co) . The phase of 

the continuous output of each radiator is randomly modulated by the dephasing effect of 
collisions which occur random at some mean collision rate of x" 1 . The figures on the 

next page shows samples of the field radiated per atom (left-hand graphs) and the 
intensity radiated per atom (right-hand graphs) for ensembles of 50, 100, 200 and 300 
collision interrupt atoms. 7 . This model is treated "experimentally' in the discussion 
entitled A "Toy" Model of a Randomly Fluctuating Optical Field which is located at 




[ VIA-21b ] 



given by 




[ VIA-22 ] 



7 



The figures show samples of the sums 



fv 1/ r v 1/ 

I^SUKP; I IN and |2^Sin(]),. | N for N atoms with the ph; 





ase 



radiated by each atom (fv distributed uniformly between and 271 . 
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http://www.deas.harvard.edu/courses/ap216/lectures/ls_3/ls3_u6A/toy/toy.html. From 
that discussion, it is reasonable to write for a collision-broadened radiator 

Y (1) (x)= exp[-i (co)x- (x/x )] [ VIA-23 ] 

(see a graph of this function below). From Equation [ VIA- 15 ], we see that the 
spectrum of a collision-broadened source is given by the Lorentzian form 

f(g>) = ^r[(<»- (®)) 2 +0K) 2 ] [ VIA - 24 ] 

Field and Intensity Fluctuations for Collision- or Doppler-Broadened Sources 



Mean =0.00494765 

Variance/ Mean =0.95 




10 20 30 40 50 
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Doppler-Broadened Sources -- An example of inhomogeneous-Broadening: 

A rather more intricate model of a chaotic source is that of a Doppler-broadening 
system of oscillators. In this model one pictures an ensemble of randomly moving 
atomic radiators each of which emits a field of constant amplitude, oscillating at a fixed 
frequency, but the observed frequency of a given atom is Doppler- shifted with respect 
to that fixed frequency. In the Doppler model, we write the sum of fields emitted by a 
ensemble of moving radiators as 

E(t) = E £exp {-/ [s>. r — 4> J} [ VIA-25 ] 

i 

where the ct), .'s are the frequencies of individual atoms which are Doppler shifted from 
the mean radiated frequency (to) . In this model, the effects of collisions are neglected 
and, consequently, the first-order correlation function is written as 8 



Since the dynamics of the individual atoms are uncorrelate, all cross-terms vanish. 
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£(V 2 ))= E 2 J £(^p{i[^{t l -z 1 /c)+^> i -co.(; 2 -z 2 /c)H).]}\ 



/ 

'• J [ VIA-26 ] 

= £?E(exp[-i'(D I . x]) 



where x= t 2 - z 2 / c- 1 { + zj c . The sum may be converted to a integral over a Gaussian 
distribution of Doppler- shifted frequencies and, hence, 



[E*( Zl t) E(z 2 t 2 ))= El J exp[-/co, x] expj-(co,. - {co)) 2 /25 2 } Jco,. 

o 

= ^V atomJ E 2 exp[-/(a))x-5 2 x 2 /2] 



[ VIA-27 ] 



The corresponding degree of first-order temporal coherence is (see a graph of this 
function below) 

Y (1) ( x) = exp [- i {©) x - 8 2 x 2 /2] [ VIA-28 ] 

and the associated normalized power spectral density follows from Equation [ VIA- 15 ] 

as 

itco) = [27i5 2 ]" ! ' 2 exr|-(co-(co))725 2 j [VIA-29] 
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Degree of First-Order Coherence - Absolute Value 




Spatial Coherence: The van Cittert-Zernike 

In the rudimentary treatment of the Young interference experiment recapitulated above, we 
assume an ideal point source and, thus, neglect any effects of spatial coherence . In 
particular, we assume that the field at points Pi and P2 due to a given radiator are in time 
synchronism. However, when we deal with extended sources, we must enlarge our notion 
of coherence The subject of spatial coherence was first developed by van Cittert and, later, 
it was extended by Zernike. 9 To define the problem, let us consider the following 
geometry: 



9 See P. H. van Cittert, Physica, 1, 201 (1934) and F. Zernike, Physica, 5, 785 (1938). 



R. Victor Jones, April 27, 2000 



The Interaction of Radiation and Matter: Quantum Theory 



Page A 15 




The task before us is to determine the so called mutual intensity J(P V P 2 ) due to the 
extended source which might be observed for the two points located on the observation 
plane. For simplicity, we assume that extended source is in a plane parallel to the 
observation plane and that it can be divided into cells of statistically independent radiators. 
Thus, the total field at any point on the observation plane is given by summing the fields 
due to each of the cells — viz. 

E(P,t)=Y,E k {P,t) . [VIA-30] 

k 

The mutual intensity is then defined as 

J{P„P 2 )={E(P v i)E*(P 2 ,t)) 

= L(£*M <(M)+LL(^(^)4(M) ' [ VIA " 31a] 

k k'*k k 
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Because of the assumed statistical independence of the cells, the cross-terms vanish so that 

J(P l ,P 2 )=L(E k {P l ,t) E*{P 2 ,t)) . [VIA-31b] 



It the spirit of the Huygens principle , we assert that component fields radiated by each 
coherent cell are spherical waves of the form 



E k {P,t) = E 



PS„ 



[expj^-/((o) t-*-^ 1 



psaV 



PS, 



!\ [VIA-32] 



where E(S k ,t) is field radiated by each element of the source and P S k is, obviously, the 
distance between the source and observation positions. Therefore, we can write 



(E k {P lt t)Er k {P 2 , t))={E 



c )/{ [PiS k ][P 2 S k ] j 



/ , , *( < x# \\ f ex p[ iM( p iS k -P 2 S k )/c] 

<E(S„,)E ( Sl ,,- ft S, - P, S,)/ C )) { \ P 2 1[PiSi] "\ 



[ VIA-33 ] 



In most instances, the path difference will be small compared to the coherence length of the 
cell, so that, in the limit of many small cells, we obtain the famous van Cittert-Zernike 
theorem in the form 



n exp[i{A)(P 1 S-i , 2 S)] 
7(5)- 



So u 



[n s] [p 2 s] 



j4 [VIA-34] 



and the first-order degree of mutual (spatial) coherence is defined as 
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, , J (PA) 
y a =y{P lt $- ' - 



** Source 



where 



Source 



2R 

and, hence, we may approximate Equation [ VIA-35a ] as 



[ VIA-35a ] 



KP n )= JJ-i%<*A, ■ [VIA-35b] 



[p„sy 

For most problems of interest, we may take 

Pfi S = J* 2 + (*„ ~ $ 2 Hy„ -Ti) 2 ~R + {X "- ^ + ( ^ - ^ [VIA-36] 



. . Source 

112 ~ 



j j [ ' ' L ^ R Jj 



[ VIA-37a ] 



Sou 



Further, we see that in the Fraunhofer or far-field approximation, the mutual 
coherence is the Fourier transform of the source intensity distribution — viz- 
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s 



l(^,T\)exp[i{k)^(x l -x 2 )/R] exp[ i(k) T\( yi - y 2 )/R] d^dy\ 



~ Wee rr ! [ VIA-37b ] 



]J/(^t|)^Jti 



Source 

— and, consequently, the measurement of the degree of mutual coherence of a 
remote source provides a means to measure the angular size of that 
extended sources -- i.e. stellar objects. 

We see then that a complete first-order description of the coherence properties of a classical 
field is embodied in the complex first-order degree of spatial-temporal 
coherence which is defined as 



y ffl {\W 2 )= n = = = [VIA-38a] 



where 



(E*(\,t)E(r 2 ,t 2 )) 

[\ E M\ 2 )\^(^)\ 2 



J 

(E^, t) E (i 2 ,t 2 )) , ± \ E%,t ] + t) E{r 2 ,t 2 + t)dt ■ [ VIA-38b ] 



Free Space Propagation of Coherence Functions: 

We derive here an equation which allows us to generalizes the notions implicit in the 
van Cittert-Zernike theorem by providing us a description of how coherence propagates 
through space. We start by assuming that the analytic signal representing the field of 
interest satisfies an wave equation of the form 
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Vj E (?„ ,t n ) = \^- 2 E (?„ ,t n ) [ VIA-39 ] 

c at,. 



If we multiply through by E (r 2 , t 2 ) , we see that 



= ^^- 2 E(i l ,t l )E%,t 2 )\. [VIA-40a] 



c 

i a 2 



Similarly 

v 2 r (?, ^ ? 2 ? 2 ) = 1 ^ r ft ? p ? 2 ? 2 ) . [ viA-40b ] 

Therefore the coherence function must satisfy the fourth-order equation 

1 3 2 

V 2 1 V 2 2 r(r l ,r 2 ;x)= r(i lt i 2 ; %) [VIA-41] 



The implication of Equations [ VIA-40 ] and [ VIA-41 ] is that coherence is a field that 
propagates through space and may be treat by methods which apply to other field 
variables. In fact, the content of the van Cittert-Zernike theorem is a description of how the 
coherence of extended source propagates to remote observation points. This point of view 
is particularly valuable in the treatment of wave propagation through random media. As a 
beam propagates through such a medium, the effect of random scattering events scrambles 
the phase of the beam and the simple notion of a field rapidly loses meaning. However, the 
evolution of the coherence function provides an important measure of the statistical 
properties of the random medium. 
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Higher-Order Correlation Functions -- A classical explication of 

THE FAMOUS HANBURY BROWN-TWISS EXPERIMENT: 

Generalizing Equation [ VIA-38a ], the degree of nth-order spatial-temporal 
coherence can be defined as 10 



r,)f) fc(i„4) (|e (?„„ qf ) (|B (£gf) 



[ VIA-42 ] 



In particular, the degree of second-order temporal coherence is defined as 



(2) — \ - - - / \ i 

Y (l * rrf ' x) ^ 77" Tx 1 = 777^ 72 [via-43] 

(Kw)| } ( 7 (w)) 

We shall see that this function is an important measure of the relative timing of of intensity 
fluctuations. 11 As a first step, let first find an expression for the average intensity radiated 
by a collection of independent oscillators — viz. 



|£(r ref ,; + x)| \ (/(r re ,r)/(r ref ,f + x) 2 ) 



10 R. J. Glauber, in Quantum Optics and Electronics, Les Houches, 1964 (edited by C. DeWitt, A. Blandin, and 
C. Cohen-Tannoudji), p63, Gordon and Breach (1965). 

11 By the famous Schwartz inequality 



</l„l 2 \/l Al 2 \ 



( f (? ref ,'K(W+^) 2 )< ( I(r ief ,t) 2 }( I (v xeV t+x) 2 ) 
and quite generally 



so that 



Y (2) (? ref ^)<l 
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/ 3 (7 (,))= [e c/2] ( |£(.)| 2 )= [E c/2] z£(f)\ 

\ u / [ VIA-44 ] 

= [e.'/2]E(|3«| 2 ) 

In this expression all of the cross-terms are presumed to vanish, since the radiation from 
any given atom is statistically uncorrelated with that of any other atom. For the collision- 
and Doppler-broadening models discussed above, the oscillators differ only in phase and 
this average can be expressed as 

i^(i{t)}=N atom [e c/2]\af. [VIA-45] 

To obtain the root-mean- square deviation in the cycle average of the intensity, we first 
calculate 

([i {•)]*)= [E, f /2] 2 ( |£(<)| 4 ) = [e,e/2] J /Ee,(<)E^KM<(<) ) 

=hciirfo\mYL{m?m\ 2 )\ 

l < i*j J 

= [e. C /2] 2 |l(|£,(')|>E(|£,(0| 2 )(|£^)| 2 )} 

I i i*j * ' J 

= |iV atom ( |£,^)| 4 ) + 2iV atom (iV atom -l) [{ ^(ttf 

Again, for the collision- and Doppler-broadening models, this average can be expressed as 



[ VIA-46 ] 



■V 2 
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([/(^)] 2 )=[^ o J 2 [e c/2] 2 |a| 4 [2-(^ tom ) 



= 2 



[ VIA-47 ] 



and, thus, we find to a very good approximation that 



P ($)-[(!(<)) 



= I 



[ VIA-48 ] 



which is consistent with the experimental results presented above. Let us now look at the 
two-time average of the intensity from a collection of statistically independent 
oscillators — viz. 



(2(02W)=[e c/2] 2 (|E(0| 2 |E(t l )r) 

\ ij.kj I 



[ VIA-49 ] 



As a result of statistical independence, most of the cross-terms vanish so that 



(7(07W)=[e c/2] 2 |j:(|B,(0| 2 [£ f fe)| 2 ) 



+ 



i( mt\ E Mt) + i{ m $m m 



[ VIA-50 ] 



1*1 



Again, for the collision- and Doppler-broadening models 
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{/(0/(^))=Kc/2] 2 {iv text (|£,(0| 2 K(rf) 



(m\ 2 ) 2+ i E M^) 



[VIA-51 ] 



and, hence, 



y» = i + |y w (t)| 2 



[ VIA-52 ] 



which is important result which holds for any chaotic source — see the figure on next page. 
Notice that this result predicts a bunching of the intensity variations. 




The Hanbury Brown - Twist Experiment 
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Schematic diagram of the Hanbury Brown-Twiss experiment 
demonstrating coherent interference in light from incoherent 
sources. Sources 1 and 2 are separate by distance d1 2 } and emit 
photons of identical wavelengths which are detected by detectors 
A and E located a distance dAB apart. The distance between the 
source plane and the detector plane L. A coincidence between 
detector outputs results in a composite signal exhibiting an 
interference effect depending on both d1 2 and dAB. 
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Light 
Source 




n 



Electronic 
Correlator 



The Hanbury Brown - Twist Spectrometer is designed to measure 12 

([Afref - A] [4(?ref.' 2 ) " h}) => "4 {('ifref ''J ^fef t 2 ))~ P) [ VIA-53 ] 



This correlation may then be written 



= y (2) (x)-l 



[ VIA-54 ] 



Generally, the results of these experiments are consistent with Equation [ VIA-52 ], but 
some of the more detailed observations require a quantum mechanical interpretation. 



12 R. Hanbury Brown and R. Q. Twiss, Proc. Roy. Soc, A 242, 300 (1957). 
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